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We investigate a (3+l)-dimensional hydrodynamic expansion of the hot and dense system created 
in head-on collisions of Pb+Pb/Au+Au at beam energies from 5 — 200 A GeV. An equation of state 
that incorporates a critical end point (CEP) in line with the lattice data is used. The necessary 
initial conditions for the hydrodynamic evolution are taken from a microscopic transport approach 
(UrQMD). We compare the properties of the initial state and the full hydrodynamical calculation 
with an isentropic expansion employing an initial state from a simple overlap model. We find 
that the specific entropy (S/A) from both initial conditions is very similar and only depends on 
the underlying equation of state. Using the chiral (hadronic) equation of state we investigate the 
expansion paths for both initial conditions. Defining a critical area around the critical point, we 
show at what beam energies one can expect to have a sizable fraction of the system close to the 
critical point. Finally, we emphasise the importance of the equation of state of strongly interacting 
matter, in the (experimental) search for the CEP. 



Heavy ion collisions at intermediate incident beam ener- 
gies (5 — 200/1 GeV) have recently attracted more and 
more attention, since one expects to be able to scan a 
wide - and by current expectations highly interesting 
[3, - region of temperature T and baryo-chemical po- 
tential [Lb in the QCD phase diagram. Following re- 
cent theoretical investigations (for recent lattice QCD 
results see [a 14 M, fo r phenomenological studies see 
d, 0, IE S fiol. Illl. 12l|). one hopes to find experimen- 
tal evidence for a phase transition from hadronic mat- 
ter (with broken chiral symmetry), to a Quark-Ghion 
Plasma (QGP) phase (where chiral symmetry is re- 
stored). Especially the so called critical end point (CEP), 
a point in the phase diagram that terminates the phase 
transition-line of the first order transition (which is ex- 
pected for high chemical potentials), is of great interest. 

Key observables like the directed and ellip tic flow v\ 

and v 2 jsiiijii, m, e s e 0, e mi e e e 

I2H I27I I28L l29l. 3C4] - but also particle multiplicities, ratios 
and their fluctuations - are of imminent interest. They 
have been predicted and sometimes already shown to be 
sensitive to the properties of the QCD matter, i.e. to the 
Equation of State (EoS) and the active degrees of free- 
dom, in the early stage of the reaction. And indeed, the 
energy dependences of various observables show anoma- 
lies at low SPS energies which might be related to the 
onset of deconfinement @, EH • 

For a hydrodynamical modelling of heavy ion reactions 
to study these observables one needs to specify initial 
conditions, i.e. an infinite set of space-time points with 
their corresponding energy- and baryon density. Since 
experimental data provides mainly information that is in- 
tegrated over the systems time evolution, the initial state 
for hydrodynamical simulations has to be inferred from 



model assumptions or by educated 'guessing' in compar- 
ison to data. The latter approach - usually applied for 
relativistic hydrodynamical simulations of nuclear colli- 
sions - is, however, by no means straight forward and 
highly non-trivial: the connections between (observed) 
final state and the inferred initial conditions is blurred 
by the unknown equation of state, potential viscosity ef- 
fects, and freeze-out problems. Another issue concerns 
the assumption of thermal equilibrium, which is proba- 
bly not true at least for the early stage of a heavy ion 
collisions at intermediate energies. 

There have been attempts to solve these problems 
by describing such collisions with viscous or multi-fluid- 
lwdrodynamic models [H, [HI, [H, HE H3, [H, [H, H3, [H, 

la ESSES!, 

aut the practical application of these 
models is difficult. To avoid (some of) these problems in 
this paper, we describe the initial stages of the collision 
with a non-equilibrium transport model (UrQMD). We 
then use the so obtained distributions for energy- and 
baryon-density as initial conditions for a one-fluid but 
fully (3+ l)-dimcnsional hydrodynamical calculation. For 
the hydrodynamical evolution an EoS with a first order 
chiral phase transition and a CEP at finite [Ib is applied. 

This paper is organised as follows: In the first part we 
describe the chiral equation of state that is used in more 
detail and explain how the structure of the phase dia- 
gram with a critical end point is modelled. Afterwards 
two different scenarios for the initial conditions and the 
subsequent evolution are compared. The full hydro evo- 
lution with microscopic initial conditions is contrasted 
with lines of constant entropy per baryon number from a 
simple overlap model. Then, the results concerning the 
CEP are presented and the influence of the equation of 
state on the evolution is studied. The last part sum- 
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marises the paper. 

Let us start with the discussion of the equation of 
state. The present chiral hadronic SU(3) Lagrangian 
incorporates the complete set of baryons from the low- 
est flavour-5'C/(3) octet, as well as the entire multiplets 
of scalar, pseudo-scalar, vector and axial-vector mesons 
44j . In mean- field approximation, the expectation values 
of the scalar fields relevant for symmetric nuclear matter 
correspond to the non-strange and strange chiral quark 
condensates, namely the a and its ss counterpart £, re- 
spectively, and further the uj and <f> vector meson fields. 
Another scalar iso-scalar field, the dilaton x> is intro- 
duced to model the QCD scale anomaly. However, if \ 
does not couple strongly to baryonic degrees of freedom 
it remains essentially "frozen" below the chiral transition 
[3] • Consequently, we focus here on the role of the quark 
condensates. 

Interactions between baryons and scalar (BM) or vec- 
tor (BV) mesons, respectively, are introduced as 



BM 



(swo- + gicO i>i 



(i) 



,7ow° + 5^7o</>°) rpi , (2) 



Here, i sums over the baryon octet (N, A, S, S). A term 
£ voc with mass terms and quartic self-interaction of the 
vector mesons is also added: 

£vec = ^a^x 2 ^ 2 + \^<$,X 2 <t> 2 + 54 + ^) • 



The scalar self-interactions are 
C = 



^oX V + C 2 ) + fc! (- 2 + C 2 ) 2 + fc 2 (y + C 4 ) 



+ k 3 x<J 2 C - k±X 4 - 7X 4 In ^4 
4 Xq 



•J CTq (,0 



Interactions between the scalar mesons induce the spon- 
taneous breaking of chiral symmetry (first line) and the 
scale breaking via the dilaton field x (last two terms). 

Non-zero current quark masses break chiral symmetry 
explicitly in QCD. In the effective Lagrangian this corre- 
sponds to terms such as 



Xo 



l l ho + (V2m 2 K f K 



V2 



(•4) 



According to £bm ©, the effective masses of the 
baryons, m*(a, () = gi a a + g^ ( , are generated through 
their coupling to the chiral condensates, which attain 
non-zero vacuum expectation values due to their self- 
interactions [44| in £ ©• The effective masses of the 



mesons are obtained as the second derivatives of the 
mesonic potential Vmcsou = — £0 — C vcc — £sb about its 
minimum. 

The baryon-vector couplings <7, w and gj& result from 



n l q)g% , g it p = —{n l s — n\)V2g% , where gj denotes the 
vector coupling of the baryon octet and n l the num- 
ber of constituent quarks of species i in a given hadron. 
The resulting relative couplings agree with additive quark 
model constraints. 

All parameters of the model discussed so far are fixed 
by either symmetry relations, hadronic vacuum observ- 
ables or nuclear matter saturation properties (for details 
see [3). In addition, the model also provides a sat- 
isfactory description of realistic (finite-size and isospin 
asymmetric) nuclei and of neutron stars 3, 3, 4(jj ]. 

If the baryonic degrees of freedom are restricted to the 
members of the lowest lying octet, the model exhibits 
a smooth decrease of the chiral condensates (crossover) 
for both high T and high /x 4J, |47| . However, addi- 



tional baryonic degrees of freedom may change this into 
a first-order phase transition in certain reg imes of the 
T-[i q plane, depending on the couplings 47, 48, 4§|. To 
model the influence of such heavy baryonic states, we add 
a single resonance with mass = mo + gaa and vec- 
tor coupling g^ — rygNu; ■ The mass parameters, mo , 
gn and the relative vector coupling ry represent free pa- 
rameters, adjusted to reproduce the phase diagram dis- 
cussed above[86j. In principle, of course, one should cou- 
ple the entire spectrum of resonances to the scalar and 
vector fields. However, to keep the number of additional 
couplings small, we effectively describe the couplings of 
all higher baryonic resonances by a single state with ad- 
justable couplings, mass and degeneracy. This method 
of replacing the influence of many states by an effective 
resonance has a long tradition in scattering theory and 
related fields (see e.g. (5pl|). 

In what follows, the meson fields arc replaced by their 
(classical) expectation values, which corresponds to ne- 
glecting quantum and thermal fluctuations. Fermions 
have to be integrated out to one-loop. The grand canon- 
ical potential can then be written as 



Cl/V — — £ VC c — Cq — £sb — V v . 



(5) 
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pure /-type coupling as discussed in [44[, g. 



where 7b, 7a/ denote the baryonic and mesonic 
spin-isospin degeneracy factors and E* B M {k) = 

k 2 + m* B M 2 are the corresponding single particle en- 
ergies. In the second line we sum over the baryon octet 
states plus the additional heavy resonance with degen- 
eracy 7b (assumed to be 16). The effective baryo- 
chemical potentials are /1* = fa — g^oj — g^cj), with 
Mi = { n \ — n q)f i q + i n l ~ n \)^s- The chemical poten- 
tials of the mesons are given by the sum of the corre- 
sponding quark and anti-quark chemical potentials. The 
vacuum energy V vac (the potential at ps = T = 0) has 
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been subtracted. 

By extremizing Q/V one obtains self-consistent gap 
equations for the meson fields. Here, globally non-strange 
matter is considered and (x s for any given T and [i q is ad- 
justed to obtain a vanishing net strangeness. The domi- 
nant "condensates" are then the a and the u> fields. There 
have also been first attempts to model the dynamical evo- 
lution of the condensates themselves instead of 'locking' 
them at their equilibrium values (see e.g. [5lLl5^|). 

Let us now turn to the explanation of the initial con- 
ditions that have been used. Two different ways to de- 
scribe the initial conditions and the expansion will be 
discussed. In one setup the energy and density distri- 
butions obtained from the UrQMD transport simulation 
are mapped to the thermodynamic quantities, which then 
serve as initial conditions for the (3+1) dimensional hy- 
drodynamic evolution. In the second setup initial energy- 
and baryon-densities as obtained from a simple overlap 
model are employed and then paths of constant entropy 
per baryon are followed. 

In the first scenario, the Ultra-relativistic Quantum 
Molecular Dynamics Model (UrQMD 1.3) is used to cal- 
culate the initial state for the hydrodynamical evolution. 
This has been done to account for the non-equilibrium 
in the very early stage of the collision. In this config- 
uration the effect of event-by-event fluctuations of the 
initial state is naturally included. Due to the early 
transition time to hydrodynamics, only initial scatter- 
ings, i.e. baryon-baryon collisions and string excita- 
tions/fragmentations are relevant here. As many details 
of the UrQMD model are not relevant for the present 
initial state calculation we refer the interested reader for 
the details 53|, 54 1 of the UrQMD model and its thermo- 



dynamic properties 55j, l56l . 157] . The coupling between 



the UrQMD initial state and the hydrodynamical evo- 
lution happens when the two Lorentz-contracted nuclei 
have passed through each other. This assures that (es- 
sentially) all initial baryon-baryon scattering have pro- 
ceeded and that the energy deposition has taken place. 
It should be noted that the present approach is different 
from using a kinetic model for the freeze-out procedure 
H [H El ED, E3, El HH which is not done in the 
present investigation, but it is in spirit similar to the 
NeXSPheRIO approach [H, EE S M, 0, [zl | • 



To allow for a consistent and numerically stable map- 
ping of the 'point like' particles from UrQMD to the 3- 
dimensional spatial-grid with a cell size of (0.2fm) 3 , each 
hadron is represented by a Gaussian with a finite width. 
I.e. each particle is described by a three-dimensional 
Gaussian distribution of its total energy-, momentum- 
fin x-, y-, and z-direction) and baryon number-density. 
The width of these Gaussians is chosen to be a = 1 fm. A 
smaller Gaussian widths leads to numerical instabilities 
(e.g. entropy production) in the further hydrodynam- 
ical evolution, while a broader width would smear out 
the initial fluctuations to a large extend. To account for 
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FIG. 1: The initial temperature distribution in the z-y plane. 
Where z is the beam-axis and y the out of plane axis. 
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FIG. 2: Excitation functions of S/A for the geometrical over- 
lap model and UrQMD initial conditions. 



the Lorentz-contraction of the nuclei in the longitudinal 
direction, a gamma-factor (in longitudinal direction) is 
included. The resulting distribution function, e.g. for 
the energy density, then reads: 

(x - x p f + (y- y v ) 2 + (j z (z - z p )) 2 



e c i{x, y, z) = Nexp - 



2a 2 



(6) 



where N — ( w-)^ -^Ei a b provides the proper normal- 
isation, e c f and E c { are the energy density and to- 
tal energy of the particle in the computational frame, 
while (x p ,y p ,Zp) is the position vector of the particle. 
Summing over all single particle distribution functions 
leads to distributions of energy-, momentum- and baryon 
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number-densities in each cell. 

This is done for Pb+Pb/Au+Au collisions at Ei a h = 
5 — 200 A GeV with impact parameter b = fin. To relate 
the distributions of energy and baryon number-density 
to thermodynamic quantities like pressure, temperature, 
chemical potential or entropy-density, the equation of 
state described above is used. As an example, Fig. Q] 
shows the initial temperature distribution obtained for 
E lah = 10A GeV. 

We contrast the microscopically calculated initial con- 
ditions described above with a simplified overlap geome- 
try initial condition. Therefore, we assume that the en- 
tire initial beam energy and baryon number equilibrates 
in a Lorentz-contracted volume determined by the over- 
lap of projectile and target in the center-of-mass frame. 
This allows to obtain a straightforward estimate for the 
initial baryon number and energy density: 




1 2 3 4 5 6 7 




2 7cms Po 



initial 

Pb 
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e = \/spo 7CMS 



(7) 
(8) 



It was shown in [35| that in the energy range of interest 
here, this rather simple approach reproduces the specific 
entropy production from a three-fluid model quite well. 

For the subsequent hydrodynamical evolution of the 
system we apply a fully (3+l)-dimensional one- fluid 
model. The hydrodynamical equations are solved by 
means of the SHASTA (SHarp And Smooth Transport 
Algorithm) as described in [z||. The EoS with a CEP is 
provided in tabulated form with a fixed step size in en- 
ergy and baryon density: Ae = 0.1 and An = 0.05 where 
e and n are given in units of nuclear ground state densi- 
ties (e w 138.5 MeV and n » 0.15fnT 3 ). This hydro- 
dynamical model has been tested vigorously and applied 
successfully for various initial conditions and physics in- 
vestigations [zi [3, 0, M, IrE S III El . 



Fig. [2] depicts the excitation function of the total en- 
tropy S per baryon number A for the initial stage of 
the hydrodynamical evolution for both initial conditions 
(UrQMD, solid circles; overlap model, dashed line). We 
have checked that both quantities are separately con- 
served throughout the whole hydrodynamical evolution, 
so S/A is a time independent constant [8 7j. Interestingly, 
the simple geometric overlap model and the UrQMD ini- 
tial conditions yield basically the same value of S/A for 
a given incident energy. 

The (isentropic) expansion paths for different beam en- 
ergies in the e- n plane are shown in Fig. [3l Here n 
is the baryon-number density. Again lines of constant 
S/A for overlap initial conditions (blue open circles and 
lines), and the (3+l)-dimensional hydrodynamical evolu- 
tion with UrQMD initial conditions are compared. The 
mean values are obtained by weighting the value of a spe- 
cific quantity in a given cell with the energy density of 
that cell. E.g., the mean number density is calculated 



FIG. 3: Isentropic expansion paths in units of ground state 
densities (eo = 138.5 MeV and no = 0.15fm~ 3 ). Red circles 
correspond to multiple UrQMD events for the same beam 
energy and impact parameter. The black region indicates the 
region below T < 0. 



from: 



< n >= 



(9) 



where i,j,k represent the cell indices. The mean val- 
ues for e, T and a q (quark chemical potential) are cal- 
culated accordingly. This has been done for equal time 
intervals of At = 2.4 frn/c. As expected, for a given 
beam energy, the dynamical paths obtained from the 
(3+l)-dimensional hydro evolution agree quite well with 
the lines of constant S/A. The variation in energy and 
baryon density due to the variation of the initial state 
in UrQMD even for a fixed impact parameter and fixed 
beam energy are studied for E\ a ^ = 10^4 GeV and indi- 
cated by the red full circles. 

As a next step the hydrodynamic evolution of the sys- 
tem is shown in the T — a q plane in Fig. [U Also in- 
dicated is the first order phase transition line and the 
CEP of the employed chiral EoS. Included are again 
lines of constant S/A and hydrodynamical evolution 
paths for the same beam energies (from left to right: 
160, 100, 40, 10, 5A GeV) as in Fig. H As one can see, the 
mean temperatures and chemical potentials of the hydro- 
dynamical evolution are not identical to the respective 
lines of constant entropy. This is due to the averaging 
procedure while a single cell does follow the isentropic 
path. The time evolution of a central cell at the origin 
in T and u q is depicted in comparison (green dashed line 
with green open circles). 

It is worthwhile to connect the present discussion of 
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FIG. 4: Isentropic expansion paths in the T — fi q plane for 
very central Pb+Pb/Au+Au reactions. UrQMD initial con- 
ditions with (3+l)-dimensional chiral hydrodynamical evolu- 
tion (averaged, full red line with circles; central cell, dashed 
green line with circles) , isentropic expansion from the overlap 
model initial conditions are shown as full line in blue. Beam 
energies are from left to right: 160, 100, 40, 10, 5A GeV. The 
phase boundary of the model is shown as full black line with 
the critical end point, the Fodor and Katz critical end point 
is shown separately with error bars [85J. 



the evolution path with UrQMD initial conditions to the 
results obtained with various models (however without 
phase transition) that are investigated in [l2| . The cur- 
rent findings with the chiral equation state support the 
main statement given there that it might be possible to 
reach the phase boundary to the QGP/chiral restoration 
already at moderate beam energies (-Ei a b ~ 5— 10 A GeV). 
A further discussion about the effect of different equa- 
tions of state will be presented below. It should also be 
noted that we do not observe a focussing of the hydro- 
dynamical trajectories towards the critical end point in 
contrast to the findings by 8l|, l82| . However, see also the 



FIG. 5: Time evolution of the critical volume for different 
beam energies for the CEP obtained with the chiral EoS. 



discussion in 8^, 84] 



Having at hand an EoS with a critical end point it is 
possible to explore, which fraction of the evolving system 
stays for how long close to the critical end point. The en- 
ergy dependence of this exposure tine is also investigated. 
Therefore, we define a 'critical volume' by adding all cells 
that have a temperature of Xcep ± 10 MeV and a chem- 
ical potential ^cep ± 10 MeV for each time step. Fig. 
shows the time evolution of this critical volume at vari- 
ous incident energies for the critical end point obtained 
by the chiral EoS. Despite the surprising fact that the 
maximal volume is reached for the highest beam energies 
(Eub = 160 — 200/1 GeV) a quite large critical volume of 
around ~ 200 fm 3 is already predicted at lower energies 
of E Uh = 60,4 GeV. 
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FIG. 6: Time evolution of the critical volume for different 
beam energies the CEP obtained by Fodor and Katz 85]. 



Fig. [5]shows the critical volume for the Tcep and /icep 
values obtained by lattice QCD calculations (85|, however 
using the chiral EoS for the dynamics. In contrast to the 
chiral values in Fig. [5] the time for the maximum does 
not change with the energy in this case. The highest val- 
ues for the critical volume are still reached at the highest 
beam energies, but one has to keep in mind that the crit- 
ical end point for the volume calculation and the critical 
point in the evolution are different. The influence of dif- 
ferent EoS on the critical volume will be discussed below. 

To pin down the beam energy to maximise the critical 
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FIG. 7: Excitation function of the critical space-time volume 
for the two different CEPs. 



volume for the longest period of time, the total 4- volume 
is obtained by integration of the critical volume over the 
time. The space-time volume is shown in Fig. [7] as a 
function of .Eiab- 

Of course the actual size of the critical volume de- 
pends on the chosen AT and A/i intervals around T c 
and \x c . Also the choice of the Gaussian width of the 
particles in the initial condition does have a small influ- 
ence on the critical volume, since at larger a the distri- 
butions for T and /j, q are wider and more 'smeared out'. 
However, we have checked that for the CEP obtained 
with the chiral EoS the excitation function is indepen- 
dent of a. The integrated critical volume (for the CEP 
obtained from the chiral model) does steadily increase 
up to E\ a b ~ 2Q0A GeV. Having the maximum critical 
space-time volume in the excitation at such an unexpect- 
edly high beam energy is due to the way the lines of 
constant entropy behave at the phase boarder (see Fig. 
EJ. With the present chiral EoS, the trajectories com- 
ing from higher temperatures turn right (towards lower 
temperatures) at the phase transition line and therefore 
much higher energies are needed as if they would turn to 
smaller \x q (left) and go up the phase line. 

As a last step, we modify our chiral EoS by adding all 
known resonances (except the baryonic decuplett, since 
it is included by means of a so called 'test'-resonance) 
with vacuum-masses up to 2 GeV as a free gas. This al- 
lows to complete the baryon resonance spectrum without 
changing the phase structure of the EoS. The position 
and type of the phase transition in the phase diagram 
stays unchanged, because the added resonances have no 
influence on the chiral condensate a. 

In Fig. [5] isentropic paths for the two chiral EoS (with 



FIG. 8: Isentropic expansion paths for different beam ener- 
gies. The dashed lines resemble a pure hadron gas EoS. The 
short-dashed and straight lines refer to the chiral EoS without 
and with the completed resonance spectrum. 

and without completed baryon resonance spectrum) are 
compared to those calculated with a free hadron-gas EoS 
(all particles have vacuum masses - which is of course not 
consistent (and only shown for discussional purposes), 
since the free hadron-gas EoS does not have any phase 
transition or CEP). As one can see, a beam energy of 
40/1 GeV is more than sufficient for an hadron gas EoS 
to reach the CEP, while the energy needed with the chiral 
EoS is 160/1 GeV (without heavy resonances) and £"i a b = 
40 — 60/1 GeV when heavy resonances are included. 



SUMMARY 

We showed, that calculating the initial conditions of 
an heavy ion collision with the UrQMD model, yields 
very similar S/A values for a given beam energy as the 
simple overlap model. We then compared an isentropic 
expansion scenario within a full (3+l)dimcnsional ideal 
hydrodynamic evolution with a chiral EoS including a 
CEP with constant S/A lines. For the hydro evolution, 
the systems mean values of energy- and baryon-density 
follow isentropic paths in the e — n phase-diagram, while 
in the T— [i plane, a single cell follows the isentropic path, 
while the averaged quantities deviate from the isentropic 
expectation. Most importantly it was shown, that con- 
cerning the search for the critical end point, it might not 
be sufficient to apply a free hadron gas EoS to estimate 
the energy needed to generate a system that, during its 
expansion, goes through the critical region. Applying 
different EoS (as we have done) can very much change 
predictions at what beam energy that CEP is reached. 
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